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Flavour changing neutral scalar interactions are a standard feature of generic multi Higgs models. 
These are constrained by mixing in the neutral meson systems. We consider situations where there 
are natural cancelations in such contributions. In particular, when the spin 0 particle has both 
scalar and pseudoscalar couplings, one may have a self-cancelation. We illustrate one such partial 
cancelation with BGL models. We also inquire whether the flavour changing quark interactions can 
lead to new production mechanisms for a neutral scalar at LHC. 
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I. INTRODUCTION 

The recent discovery of a scalar particle by ATLAS p] and CMS |2j leaves a fundamental question unanswered: 
given that a fundamental scalar exists, how many fundamental scalars are there in Nature? Indeed, within an 
SU(2)l x U(1)y gauge theory, the number of gauge bosons is fixed. In addition, measurements of the invisible width 
of the Z boson at LEP |5] have fixed the number of fermion families to be three (below Mz/ 2). So, only the number 
of spin zero particles remains to be determined. 

The most natural simple extension of the standard model (SM) is the two Higgs doublet model (2HDM) for 
a review, see for example 00. This model is interesting in itself, but also as a toy model for a wide number of 
features that could appear in other more complicated settings: among others, it includes the need to restrict the 
parameter space, such that the vacuum does not violate charge; the appearance of extra scalar and pseudoscalar 
particles; the possibility that CP is violated in the scalar sector, either explicitly or spontaneously; the existence of 
charged scalars which, like the W ± , would change flavour; and the possibility that there are flavour changing neutral 
scalar interactions (FCNSI). In this article we explore the last feature, considering two consequences. If there exists 
a scalar with FCNSI, this will have a twofold effect. On the one hand, it could conceivably be produced at the LHC 
through quark level interactions such as uc ., ds , ds, etc... On the other hand, these production mechanisms must be 
constrained by the fact that the same couplings could originate FCNSI in neutral meson systems, such as K° — K 
or B° — B , which we denote generically by P° — P . 

Flavour changing neutral interactions, such as K° — K and K° have played a crucial role in the history 

of Physics. For example, CP violation was first seen in the K° — K system [6], while the charm quark was invented 
to curtail a large contribution for K° —i R + from the box diagram with the up quark [7j. More recently, the first 
evidence of CP violation outside of the kaon system was found by the Babar and Belle experiments 00, showing 
that it arises from a large (not small) parameter, and the first evidence of Bd, s —> R + p~ has been found [TP] , 

Bounds on P° — P° mixing lead to constraints on the coupling cs q2qi of the scalar S with quarks q± and q 2 that 
scale linearly with the mass ms of S: 


CSq2qi ^ 'S’ 

Limits also arise from bounds such as B s —> fj> + /i _ , scaling quadratically as 

1 ™2 

< ' 5 '(/ 29 i m S■ 
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Generically, Eq. ([2]) provides a looser bound than Eq. (JT]) for a sufficiently large ms- Moreover, Eq. ([2]) involves the 
coupling csnn of S with p + p~, which might even be zero if the scalar S has no overlap with the (would-be) SM Higgs 
field H sm . In both cases, scalars with larger masses have a looser constraint on cs q2qi and, thus, could conceivably be 
produced at LHC via pp —> <7291 —> S. This analysis points to a (further) interesting complementary feature between 
Flavour Physics, such as would be pursued at a Super B-factory, and the scalar search to be continued at LHC’s 
Run2. Of fundamental interest is the question: how would pp —> < 72(71 —> S compare with the glue-glue production 
mechanism pp —► gg —> SI 

In Section [TTJ we introduce the mechanism of self-cancelation, possible when there are both scalar and pseudoscalar 
couplings. This is then discussed in increasingly particular cases of the two most general 2HDM, the CP conserving 
2HDM, and the BGL models. In Section |Hl| we turn to the possibility that the scalar is produced at the LHC through 
FCNSI couplings, and we conclude in section [TV] For completeness, the appendix contains formulas for FCNSI effects 
obtained for the most general Lagrangian, other than those derived in the main text. 


II. SELF-CANCELATION IN NEUTRAL MESON MIXING 
A. Generic scalar contribution 


Let us consider a spin 0 particle S interacting with two quarks qi and Qj according to 

- C Y = S {<77 (dji + ibji'js) qt + q t (a| 4 + ibZ. 75) qj j -|- ( 3 ) 

No sum on i and j is implied. If qj = < 7 ^, then a and b are real; otherwise, they are complex. In writing Eq. (|3| we 
have already used hermiticity in the form 


a ij = ,l *r kj = b-ji- (4) 

If all bij (atj) were zero, then S would be a pure scalar (pseudoscalar). Otherwise, S will be a mixture of scalar and 
pseudoscalar, and there is P violation. It is interesting to note that there are still viable models in which the 125 GeV 
scalar found at LHC has a pure scalar couplings to the up quarks, while it has a pure pseudoscalar couplings to the 
down quarks HU. 

We are interested in the neutral meson systems constituted by P° ~ qjqi and P° ~ qiqj- Cy contributes to an 
effective Hamiltonian, mediating the mixing 1121 

M 2 i = 


We denote by mp, Amp, and fp the average mass, the mass difference, and the form factor of the P Q — P° system, 
respectively. Under reasonable approximations [12j . 


(p° \n oS \ p° 

fp TOp 


24 m| l 






(5) 


A m,p = 2|Af 2 i|- 


( 6 ) 


In the vacuum insertion approximation, discussed in detail in appendix C of Ref. m , we find 


K :T = 1 


K b = 1 


m. 


('m , 

11 


(m qj + m qi ) 2 


(7) 


for the scalar and pseudoscalar operators, respectively. 

Although our main points do not depend on the exact values of and Kj *, we show in Table |I] a rough estimate 
based on the vacuum insertion approximation of Eq. |7| and the relevant input parameters. The quark masses are 


1 Notice that there are a few sign misprints in the hardcover edition of m, corrected both in the paperback edition and here. 
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Meson system 

O 

1 

O 

dO dO 

^dT^d 

dO dO 

S s 

D°-D° 

V 

1.0 

0.1 

0.1 

1.0 

Amp 

3.484 x 10 -15 

3.337 x 10" 13 

1.17 x 10 -11 

6.58 x 10~ 15 

mp 

0.497 

5.28 

5.37 

1.86 

fp 

0.156 

0.190 

0.225 

0.232 

m qi , m q2 

m d ,m 3 

m d , m b 

m s ,m b 

m u ,m c 

Ka 

25.8 

2.58 

2.56 

3.12 

Kb 

273.8 

18.3 

18.2 

24.3 


TABLE I: Meson mixing input; Amp, Mp, fp, m q in GeV. Taken from Ref. [Hi- 

taken (in GeV) as m u = 0.0023, = 0.0048, m s = 0.095, m c = 1.275, and mb = 4.2 [17. We notice that the 

pseudoscalar matrix elements K b are always larger than their scalar counterparts K%. For example, in the B systems, 
we get K b /Kp ~ 7 in the vacuum insertion approximation. It is interesting to compare with the values obtained 
from lattice. For instances, we find 



5/iXi + 6 /i B S. 


( 8 ) 


where the right-hand side involves the quantities introduced in Ref. mi 2 . In particular, B^ (q = d, s) is the bag 
parameter common to the operators O 2 = (bj^q ) 2 and O 2 = (b^Rq) 2 , with jl,r = (1 T 7s)/ 2, while B^ (q = d,s) 
is the bag parameter of the operator O 4 = (b r YLq)(bjRq). The results obtained from lattice differ from those in the 
vacuum insertion approximation by at most a factor of three. Of crucial importance is the ratio K b q /K a q ~ 3 
obtained from lattice. The fact that the ratio obtained using the lattice results of [18] is closer to unity than that 
obtained in the vacuum insertion approximation will be of interest in the following. 

Next we highlight one of the main points in this work. In Eq. |5]) the scalar and pseudoscalar components appear 
independently; there are terms in a-.. and b 2 ., but no term in One would say that they do not interfere. 

But, because aP and bf are complex, Eq. (|5^ shows the very interesting feature that the two terms can cancel each 
other. And, as we will illustrate below, there are generic classes of models in which they could easily arise with the 
opposite sign. This has the result that the scalar contribution of a spin zero particle may cancel the pseudoscalar 
contribution of that same spin zero particle. There are known instances where the contribution of some scalar cancels 
the contribution of some other scalar. This occurs, for instances, in the scalar contributions to the electric dipole 
moment of the electron in CP violating two Higgs doublet models, when in the decoupling limit m- But the feature 
present in Eq. (|5j) is something else. It is the two components of a single scalar that may cancel each other. We denote 
this effect by “self-cancelation”. As far as we know, this feature hasn’t been properly appreciated before, especially 
in the context of its implications for the LHC. 


B. Self-cancelation in the most general two Higgs doublet model 

The Yukawa interactions of quarks in the two Higgs doublet model (2HDM) may be written as 


— Cy = <Zl (Ti^i + r 2 <!> 2 )?Ti? + 9 l(A il>i + A 2 I> 2 )pr + h.c., (9) 

where 4> a (a = 1,2) are the Higgs doublets, qL = (pL,nL) T is a vector in the 3-dinrensional family space of left- 
handed doublets, and ur and pr are 3-dimensional vectors in the right-handed spaces of charge —1/3 and +2/3 
quarks, respectively. The complex 3x3 matrices Ti, T 2 , Ai, and A 2 contain the Yukawa couplings. 

After spontaneous symmetry breaking, the fields acquire the vacuum expectation values (vevs) v\/y2 and c 2 /v2- 
In general, rq and r; 2 are complex, but one may, without loss of generality, choose a basis where iq is real and 


2 Notice that our conventions for matrix elements differs from those in Ref. m by a minus sign and by l/(2mg). Of course, physical 
results are the same and, when all is properly taken into account, one obtains Eqs. § 








4 


V2 = 


= \v 2 \e l5 ■ It is convenient to perform a unitary transformation into the Higgs basis {.Hi,# 2 } through [lfjj 


H 2 


= U t 



where 


( 10 ) 


tT = - diag (l, e lx ) 
v 



( 11 ) 


is unitary, and v = \/v\ + v\ = ( \f2Gp ) 1 / 2 . As is clear from Eq. (ITT]) , all the vev is now in H\. And, since H 2 has 
no vev, its phase can be altered at will. We will choose X = 0, Ref. [21] has x = — 5 + 7 r. We may write 


G+ 

(v + H sin + iG 0 )/V 2 ’ 


H 2 


H+ 

(R + iI)/y/2 ’ 


( 12 ) 


where G + and G° are the would-be Goldstone bosons and H + is the charged scalar. H sm , R and I are neutral fields. 
The Yukawa couplings in the Higgs basis are given by = TbUba, A^ = AbU£ a . Since only Hi has a vev, the quark 
masses arise solely from Tf^ and Aj 7 as 


v r/t r H TT, 

1 Ud R 

_T_ r/t a h tj 


D d = diag (m d , m s , m b ), 
D u = diag (m u , m c , m t ), 


(13) 


where the transformations U a with a = cLl, dn,UL,UR bring the quarks into their mass basis. The couplings to H 2 
become 


^flT»U dR = N d , 

■^Ut L A»U UR = N u , (14) 


and the Yukawa lagrangian becomes 

- -j=Cy = (■ u L V, dL^DaHi + N d H 2 ) d R + (u L , J L H t )(^„H ' 1 + N U H 2 ) u R + h.c. , (15) 

where V = Ul L U dL is the CKM matrix. In general, N d and N u are not diagonal and, thus, responsible for the FCNSI 
involving R and I. Notice that H sm couples to the quarks proportionally to their masses, as would the SM Higgs 
particle. But, in general, H sm is not a mass eigenstate. 

If ViV 2 7 ^ 0, H sm and R are guaranteed not to be mass eigenstates and, if there is CP violation in the pure scalar 
sector, then I is also not a mass eigenstate. They mix through 



= T 



(16) 


into the S/ c (k = 1,2,3) mass eigenstates corresponding to the three neutral spin 0 particles. Using Eqs. (15) and 
(16), we can finally write the Yukawa lagrangian of the neutral scalar fields as 


- C Y = Y, S * { d i Ad ’ k + *75 B u ’ k ] d + u [A u ’ k +i 75 B u ’ k ] u} , 


k =1 


(17) 
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where 


and 


vA d ’ k = T lk D d + T 2k X d + + iT 3k X d _, 
vB d ’ k = -iT 2 k X d _+T 3k X d , 
vA u,k = TlkDu + T 2 k xi-iT 3k A'“, 
vB u ’ k = ~iT 2 k Xl~T 3k Xl, 


. a _ N a ±Nj 
± 2 


( 18 ) 


(19) 


for a = u,d. The down type couplings agree with those found in Eqs. (22.73)~(22.74) of [12] ■ Notice that the matrices 
A a - k and B a ’ k are hermitian, as needed for a hermitian lagrangian. 

We may now calculate the contribution to the P°-p° mixing matrix element M 2 \ in Eq. ([5| as 


24 

fpmp 


M k 


1 

m s k 



where 


( 20 ) 


»4 = 

( 21 ) 


and a = d for the K , B d , and B s systems, while a = u for the D system. 


C. Self-cancelation in a CP conserving pure scalar sector 


Although it may seem counterintuitive, one can have a spin 0 state which arises out of a CP conserving Higgs 
potential, but which, nevertheless, couples with quarks through both scalar and pseudoscalar components. To illustrate 
this mechanism, let us consider a model which is CP conserving in its pure scalar sector. By this we mean that all 
couplings in the Higgs potential are real and that both vevs are real. One may write 


Vi = VCp, 

v 2 = vsp, (22) 


implying that tp = v 2 /vi, where, thenceforth, eg, sg, and tg represent the cosine, the sine, and the tangent of some 
given angle 8 , respectively. In such cases, and continuing to consider only the scalar sector, there is one CP odd state 
(A = I), and it is common to define the lighter (h) and heavier ( H) CP even states by [221 




(23) 


where Re$° = (p a + v a )/y/2 (a = 1,2), 
this case, Eq. (161 becomes 


and in going to the second line, we have used Eqs. 



and ( |22] ). In 


( 24 ) 
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Thus, the scalar and pseudoscalar couplings of each mass eigenstate to the down type quarks become: 

V A d,H = C a -/3 Dd + S a -/3 Xf 

v B d,H = —is a -pX d , 
v A d,h = -s a -p D d + c a -p X d , 
vB d ' h = -ic a _pX d , 
vA d ’ A = iX d , 
vB d ’ A = X d . 

Similarly, for the up type quarks we find 


vA u ’ H 

— Ca—/3 D u T S a —fd 

\rU 

A + ’ 

vB u ' h 

= —i s a -/3 A'“, 


v A u ’ h 

— Sa—/3 D u T C a — 

\rU 

-P x +i 

vB u ’ h 

= -ic a _pX™, 


vA u ’ A 

= -iX", 


vB u ’ A 

= -X“. 



(25) 


(26) 


Let us concentrate on the neutral meson systems with down-type quarks. Since Dd is diagonal, the relevant 
coefficients are simplified into those listed in Table [il] Let us assume that the corresponding elements in X d and X d 


scalar 

H 

li 

A 


(xii 

- 


m 


-sl-e 

-cl-, 


TABLE II: Scalar and pseudoscalar couplings present in 2HDMs with CP conservation in the Higgs potential and in the vevs. 


have the same phase. As we will illustrate in Section [II E[ this is a rather common feature. In that case, the opposite 
signs appearing in the two columns of Table [II] imply a cancelation. Indeed, for each scalar particle (for each row 

in Table II), the contribution has the opposite sign to the (bjfj contribution. This means that the scalar 

contribution of one spin 0 particle tends to cancel the pseudoscalar contribution of the same spin 0 particle 3 . 

In getting to Table |TT| nothing was assumed besides CP conservation in the Higgs potential and in the vevs. So, the 
self-cancelation is a generic feature of these 2HDMs. For the self-cancelation to be complete in Eq. (21) one would 
need 


for h and H , while 


(K, 


-« 


(27) 

( x - 

)% k - 

-(*+: 

l), AT 

(28) 


would be needed for A. If all masses were of the same order, then a cancelation in H and h would imply a non¬ 
cancelation in A, and vice-versa. 

But there are other possibilities. Recall that s a -p controls the coupling of h to the vector bosons ZZ and W + W~. 
If h coincides with the scalar found a LHC, then s Q _ / 3 should not differ much from unity. In that case, the (? a _a factors 
in Table |TT| curtail the h contributions. Then, one could have a self-cancelation in H and a small A contribution due 
to Til a rriH, or vice-versa. 


3 


Notice that this is completely unrelated to any further cancelation which might occur between different spin 0 particles. For example, 
if one takes m h = m^ = m a and = (X^)^, then the H , h , and A contributions cancel exactly. 
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D. A conundrum of scalar-pseudoscalar mixing with a CP conserving Higgs potential 


In Section IIC we considered models where there is CP conservation in the Higgs potential, which remains unbroken 
by the vevs. In such cases, at tree level, the spin 0 states are eigenstates of CP defined in the pure scalar sector: H 
and h are CP even, while A is CP odd. Nevertheless, each spin 0 particle couples to quarks as in Eq. ([3]), meaning 
that it has both scalar (a 3 i) and pseudoscalar (bjj) couplings to quarks. Is there a contradiction? No! There is no 
contradiction. 

Let us start by considering the diagonal couplings. The point is that there is no CP violation in the pure scalar 
sector. This means that the parameters of the Higgs potential are real and so are the vevs. As a result, the tree level 
mass matrix for the neutral scalars is block diagonal and there is no CP violation in the pure scalar sector. This 
can be seen in a basis independent fashion through the basis invariant measures of CPV introduced by Lavoura and 
Silva TfjJ. They all vanish. So, where does the CP violation in ai,b lt ^ 0 come from? As is obvious from Eq. ([3]), it 
comes from the couplings with quarks; from the Nj, and N u matrices in Eqs. (14), originating in the complex Yukawa 
matrices and driving the FCNSI. Botella and Silva m have developed basis invariant measures of CP violation which 
measure precisely this type of CP violation arising from the beating of the scalar sector against the Yukawa sector. 
And the relation between these invariants and Eq. is discussed in sections 22.9.2-22.10 of Ref. [12] , What does not 
seem to have been appreciate then is that such effect can lead to self-cancelations, thus hiding potentially interesting 
FCNSI phenomena. In this respect, we stress the results found so far. The assumption that there are no cancelations 
is far from natural. It turns out that reasonable models lead naturally to (at least some degree of) cancelations. 

Let us now look at non diagonal couplings. Once in the scalar and quark mass basis, the most general CP 
transformations can be written as 


(CV) Ql (CV)' = e«* 7 °CqT, 

(CV) qj (CV)' = 

(CV) S (CV)' = r) S S, (29) 

where £ are spurious phases brought about by the CP transformation m, and we have considered only signs r]s = ±1 
for the scalar field. These can be combined into 

(CV) Sq jQi (CV)' = ns< AC ' :) Sq, qj 

(CV) Sq 0 75 ft (CV)' = -VseK'-MSqrfqj, (30) 

and Eq. ([3| is transformed into 

- (CV) C Y (CV)' = rjsS (aji - %;7s) Qj + < <l.i Of* - 7 s) 9*} d- (31) 

For CP conservation to hold, the first term of Eq. ([3| must equal the second term of Eq. (31), leading to 

a~ji = qs<i*,< l ‘ X ' £<) , 

b 3l = -vsb (32) 


Of course, the spurious phase (£ ? — £*) can be chosen to make either equation hold. This is a reflection of the known 

fact that a term by itself cannot lead to CP violation; one needs always the beating of two terms. However, these 

equations taken together mean that CP conservation implies 

= ( 33 ) 

i.e. Re [a^b-.^j = 0, which does not depend on the spurious phases. This is a rephasing independent sign of CP 
conservation. Conversely, 

Re ( a jibji) 7 ^ 0 => CP Violation (i^j). (34) 

Notice the curious possibility that one could have CP conservation with ^ 0 as long as the two couplings were 
relatively imaginary. We know of no model for which this is a compulsory feature, but the possibility should be kept 
in mind. Of course, since hermiticity of the Lagrangian requires the diagonal couplings to be real, 




CP Violation (diagonal). 


(35) 





A similar analysis for the parity transformation would lead to 


a jib]i 7^ 0 => P Violation, 

regardless of i = j ot i ^ j. 


(36) 


E. Self-cancelation in BGL models 

We have mentioned that there are models where the corresponding matrix elements of Xf and Vf have a common 
phase. One such example is provided by a model proposed by Branco, Grimus, and Lavoura, known as the BGL 
model [211 Elj- The model was constructed to obviate constraints on FCNSI by relating the matrices N u or Nd with 
off-diagonal CKM matrix elements, which are known to be small. As shown in Ref. [25] under some assumptions, 
BGL models provide the only possible implementation in 2HDMs of a relation between FCNSI and the CKM matrix 
which uses abelian symmetries. There are six such models in the quark sector 4 . Three models, known as up models 
(types u, c, and t), have a diagonal N u and a non diagonal N d - Three models, known as down models (types d, s, 
and b), have a diagonal Nd and a non diagonal N u . 


1. Up models 

After some calculations, we find for the type t model 


' m u tp 0 0 

iV„ = | 0 m c tp 0 

0 0 — m t t~] 


X u — N 

^- 1 - — ly Ul 


X u = 0, 


(37) 


( 


Nd = 


m d 


(l-\V td \ 2 )tp-\ v td |V 

—m d {tp + tg )Vt d V ts 
-m d (tp +tp 1 )V td V t l 


m. 


—m s {tp + tp l )V ts Vt d 

~(l-\V ts \ 2 )tp-\V ts \% 
-m s (tp +tp 1 )V ts V* b nib 


-m b {tp +tp l )V tb V* d ^ 
— m b (tp + tp )Vt b V ts 
(l-\V tb \ 2 )tp-\V tb \Hp 1 ] J 


XI 


( m d [(1 - 1 ^ 1 % - \ Vtd \% 1 ] - rr ^( te + tp 1 ) V ts V* d 

+ t -i) Vtd V* s m s [(1 - \V ts \ 2 )tp - \V ts \H~ x 


- rr ^(t f} + tp 1 )V tb V t * d 

+t?) Vtb VZ 


\ 


Z, ' r- / -- v / t pj £ \ * i/tj 

K + t-^)V td V? h + tp^VuVZ m b [(1 - |V th | 2 )^ - |V t6 | V 


x_ = 


m s —m d , 


2 '"{tp+tp )VtdV t % 

I! ^(tp + tp 1 )V td V t l 


2 '-(.tp + tp )VtsV t *d - 

0 


m b —m d 

2 

m b — m s 
2 


m b —m s 

2 


(tp + t^)V ts V* 


(tp + tp )V tb V t y 

(tp + tp^VtbV* 

0 


(38) 


(39) 


(40) 


The matrices A“ are hernritian, while X° are anti-hermitian. Moreover, as announced, (X^j. has the same phase 
as (pcL)j i - Notice that the (tp + tp 1 ) prefactor in the off diagonal terms can make the FCNSI large for sufficiently 
large values of tan /?. 

The matrix N d is highly hierarchical, both due to the masses and due to the CKM matrix elements. Concentrating 
on the off-diagonal matrix elements, and taking out the (tp +tp 1 ) prefactor, the orders of magnitude for X+ and X 4 
in the t type model are 


X± 


x m s 
x x 


l x x 


A 5 TOf,A 3 \ 

m b X 2 (t type), 

x / 


(41) 


4 These branch into more possibilities once one takes the leptonic sector into account m- 
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where A is the expansion parameter in the Wolfenstein parametrization of the CKM matrix m- We show only the 
12, 13, and 23 elements because we are focusing on flavour violating transitions and because the transposed elements 
are of the same order of magnitude. Thus, in the t type model, the largest contribution would occur in the B s system. 
The only difference in the c type model is that the CKM combinations V a iV*^ with a = t get changed into V a iV*A 


X ° 4 


A'| 


Thus, 


A 

m s A 

m b A 3 

y 

X 

m b A 2 

V 

X 

X 

(x 

m s \ 

m b A 3 

X 

X 

m b A 4 


X 

X 


(u type). 


(42) 


(43) 


Eqs. (41)-(43) can be used by model builders to increase some FCNSI of interest and suppress others. 


Let us focus again on the t type model. In this model 


xl 

xl 


ll 0 Q 2 ' °Qi 
TYl ao TTT-ai 


> l. 


(44) 


Since we see from Table □that K b > K q , the cancelation in Eq. (28) is not possible, while, due to the hierarchical 


mass structure of the different quark families, the cancelation in Eq. (27) is only partial, as we see from Table III 
Thus, although the relative minus signs in Table [Tl] indicate some self-cancelation, in the BGL models this cancelation 


Meson system 

O 

1 

o 

* 

dO nO 

rO pO 
*3 s ^3 s 

K?/K? 

10.6 

7.1 [2.9] 

7.1 [3.2] 

{m q2 + m qi ) 2 / (m q2 - m qi ) 2 

1.2 

1.0 

1.1 


TABLE III: Comparison of the ratios needed for the self-cancelation in Eq. (271. The numbers are obtained from Table [I] 
except those whithin [ ], which are obtained from Eq. © and the lattice simulation of Ref. un¬ 


is not complete because the ratios os masses are not enough to offset the ratios of hadronic matrix elements. We note 
that the cancelation is more effective when the matrix elements are estimated with the lattice results of Ref. HU, 
then when they are estimated in the vacuum insertion approximation. We conclude that this mechanism must be 
taken into account in the experimental search for the features of a generic 2HDM. Indeed, as this example shows, an 
accidental cancelation is quite likely and should not be ruled out a priori. 


2. Down models 


After some calculations, we find for the type b model 

*m u (1 - \V ub \ 2 )tp ~ \V ub \ 2 t^ n 

/ . 1 \ “w 


N u = 


V 


~rn u (tp +tp 1 )V ub V* b 
- m u {tp + t^)V ub V* b 


m c 


—m c {tp +tp 1 )V cb V* b 
\l-\V cb \ 2 )t 0 -\V cb \H/ 
-m c {tp + tp 1 )V cb V* b 


X“ = 


^m u (1 - \V ub \ 2 )tp - \Vub\Hp 1 


V 


- m ^(t p + t /) Vub v: b 

-^^{tp + t-^) VubV ; b 


■* 

ub 


m c 


- 0 _^ {tp+t -l )VcbVx 

- \Vcb\ 2 )t p - \ V cb \ 2 tp 
- r ^(t p + t/) Vcb V t * b 


( 0 _rn ^ {t0 + t -l )Vcb y :b 

X- = “^(t/3 + t/)v ub v* b 0 

\' r ^( t p + t/)v ub v t i + 


-m t [tp + t^^th V* b ^ 

-mtitp +tp 1 )V tb V* b 
mt (1 - \V tb \ 2 )tp - IVtb^tp 1 j 

- m ^(t 0 +t- 0 l )v tb v: b \ 

-m±rru {tf) + t- l )V tbV ^ 


mt 


a-Wto-itkiv 


mt—m 

2 

mt—vn c 


(tp + t^)V tb V^ 
f^(t p + t/)v tb v; b 
o 


(45) 


(46) 


(47) 
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N d = 


^m d tp 0 0 

0 rn s t f 3 0 



X% = N d , X d _ = 0. 


(48) 


0 0 —rribtp 

In terms of orders of magnitude, we have 

x m c A 5 mt A 3 \ 


m t X 2 J (b type), 


(49) 


xa 


x m r 

X X 
L X X 


A m t A 3 \ 
m t A 2 


(s type), 


(50) 


X\ 


^x m c A mt.X 3 ^ 
x x to* A 4 
yx x x 


(d type). 


(51) 


Notice that the orders of magnitude of the CKM coefficients in Eqs. (49)-(51) reproduce those in Eqs. (41)-(43), 


respectively; the masses are, obviously, different. Again, this is not enough to make the self-cancelation fully effective 
for BGL models in the D system. 


III. FCNSI-INDUCED HIGGS PRODUCTION AT THE LHC 

The Lagrangian in Eq. ([3]) induces pp —» S direct production with qj and qi as partons (and also, with partons q, 
and qj). At leading order in the narrow width approximation, we find 

<* \PP(9jH) ->■ S)] = 2 gjV ^ 2 (\a-ji \ 2 + \b 3i \ 2 ) Iq jQz , (52) 


where 


L QjQi 


= [ d,xf? j (x,Q 2 )f%.(x 0 /x,Q 2 )-. 

Jxn X 


(53) 


A factor of two was explicitly included in Eq. ( |52| ) (due to the fact that a given parton can come with equal probability 
from either proton - things would be different in a pp collision), N c is the number of colours, E is the energy of the 
colliding proton, /y are the relevant parton distribution functions (PDFs), and 


Xo = 


l/ln 

4E2' 


(54) 


Notice that 


g [pp(<liQj) -t S)] = Iq iqj 
° \pP(QjQi) -A S )] Iq jqi 


(55) 


is not unity, since, in general /|. ^ /£. and ^ /£. For instances, although Eq. 0 implies that \a Su \ = |«ucl and 
\b 5u \ = \buc\ > cr [pp(cu) —> 5)] does not equal a [pp(uc) — > 5)]. 

In order to estimate the change obtained in going from LO to NNLO, we use SusHi (version 1.5.0) {2E| with the 
factorization and running scales pf = Pr = ms/4 [291436] in order to find 


R{m s ) 


a [ppjbb) -> S)] n ™ L q 

(7 [pp(bb) -A S )] LO 


(56) 
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y/s = 8 TeV 


y/s = 14 TeV 




FIG. 1: Cross section for S production in pp collisions at 8TeV-Left panel (14TeV-Right panel) through the qiqj partons 
indicated, assuming that they couple as the SM bb. For reference, we show also the cross-section for glue-glue production. 


For each scalar mass ms, we use R(ms) as a universal rescaling factor for all our production cross sections, thus 
taking into account the major factors appearing in going from LO to NNLO. To see the effect of the different PDFs, 
we show in Fig. [l] the cross sections calculated as a function of ms, assuming that all couplings coincide with the SM 
bb couplings a bb = m b /v and b bb = 0. Also shown is the SM gluon-gluon fusion pr oduction cross section, typically 
two to three orders of magnitude larger than the SM bb. Of course, as seen in Fig. [lj if the FCNSI couplings were all 
equal to m b /v, then the fact that the u and d PDFs in the proton are larger than all others, would mean that the 
productions through uc, ds, and db would be the largest. We see also that there is a very small difference between sb 
and bs. 

Since the bounds from leptonic decays depend also on the details of the leptonic sector, we will concentrate on the 
bounds from P° — P° mixing. Let us decide that the scalar contribution to A/ 21 I in Eq. ([5| is a fraction 77 of the total 
contribution, including also the SM box diagram: 


MfJ = 77 \M 12 \. 


(57) 


In the notation of the CKMfitter group US CEE 

+ Mf 2 = Mi 2 = M[% A. 


Thus 


77 = 


A- 1 
A 


(58) 


(59) 


In the notation of the UTfit collaboration, A = Ce nu. For the B systems, a very conservative guess would be 
10%. In the K and D systems, the long distance contributions have a large uncertainty which could easily hide a 
scalar contributions amounting to 200% of the SM contributions, but with the opposite sign. 


A. Pure scalar or pseudoscalar couplings 


We start by assuming that the scalar S is either a pure scalar (b = 0) or a pure pseudoscalar (a = 0). From Eqs. 
([ 6 ]), and (57), we find the upper limit 


c “ = 


/ 12 11 Amp ms 
K? m P f P ’ 


(60) 


where c = a,b, and which depends linearly on ms, as announced in Eq. ([Tj). 

We have mentioned in connection with Table [I] and Eq. ([ 7 ]) that the pseudoscalar matrix elements K b are always 
larger than their scalar counterparts (at least in the vacuum insertion approximation and the lattice estimate of 
Ref. [T 8 ] used here). This means that the maximum allowed values for the pseudoscalar couplings ( b ) will always be 
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Meson system 

O 

1 

o 

* 

nO dO 

Hdr^d 

dO nO 
s ■t's 

D°-D a 

|a max | (125GeV) 

4.58 x 10" 5 

1.13 x 10“ 4 

5.61 x 10 -4 

6.28 x 10 -5 

|6 max | (125GeV) 

1.41 x 10 -5 

4.24 x 10“ 5 

2.11 x 10" 4 

2.25 x 10 -5 

|a max | (ITeV) 

3.66 x 1CT 4 

9.02 x 10“ 4 

4.49 x 10 -3 

5.03 x 10" 4 

|6 max | (ITeV) 

1.12 x 10" 4 

3.39 x 10“ 4 

1.68 x 10 -3 

1.80 x 10 -4 


TABLE IV: Meson mixing constraints on couplings for ms = 125 GeV and for ms = 1 TeV. For comparison, the bbS coupling 
is ^ = 1.138 x 10~ 2 , for a running mass of mb = 2.8 GeV. 


smaller than the corresponding scalar couplings (a) by a factor of roughly 3 (y/3 using Ref. JE]). Of course, since 
we are using estimates of these matrix elements, all results must be taken as in dica tive rather than tight constraints. 
For to/ j = 125GeV (mh = ITeV), we find the maximum values shown in Table [rv| For comparison, the bbh coupling 
is = — = 1.138 x 10 -2 , for a running mass of to& = 2.8 GeV. The best case occurs in the B° s -B° s system with 
scalar coupling, but still 


„max| |„max| 

-*^=0.05 (125GeV), = 0.4 (ITeV). (61) 

C Sbb C Sbb 


Fig- i shows the cross sections for production through q\q 2 and q 2 q\ obtained when the couplings are pure scalar 
(b = 0), having the largest possible magnitude consistent with meson mixings. Recall that the value for |a max | 
increases with ms- Similarly, Fig.[3]shows the cross sections for production through q\q 2 and q 2 q\ obtained when the 

\/s = 8 TeV -/s = 14 TeV 




FIG. 2: Sum of the cross sections a(pp (qiq 2 ) —> S) + a(pp (q 2 qi) — > S) for the maximal values of the pure scalar couplings 
consistent with meson mixings, for 8TeV-Left panel (14TeV-Right panel). 

couplings are pure pseudoscalar (a = 0), having the largest possible magnitude consistent with meson mixings. Now 
one needs the value for |& max |, which also increases with ms- 

We stress that our concern here is not on exact values. Our values for i] and the hadronic matrix elements are 
taken just to illustrate the order of the effects. As we will see in the nest section, the mechanism of self-cancelation 
can enhance the allowed effects and, thus, it should be taken as a possibility in phenomenological searches. 


B. Self-cancelations and fine tuning 

We now inquire how much fine tuning might one need in the self-cancelation, in order that a [pp(qjqi) —> S')] might 
be a relevant portion of the production mechanism. For definiteness, we concentrate on the B s system. For each value 
of ms, we define 


a \pp(bs + sb) —> S)] 


NNLO 


R a [pp(bs + sb) -A S)] 


LO 


^ [pp(gg) 

TTR , _ 1 2 

4 N7EP l|a5sl 


5)] 


SM 


\k 


L bs 


° [pp(gg) 

f Isb 


^)] 


SM 


bs I 


a- [pp(gg) S)] 


SM 


(62) 
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ifs = 8 TeV y/s = 14 TeV 




FIG. 3: Sum of the cross sections a(pp (qiQ 2 ) —> S) + a{pp (q 2 qi) —> S) for the maximal values of the pure pseudoscalar 
couplings consistent with meson mixings, for 8TeV-Left panel (14TeV-Right panel). 


In many circumstances a [pp(gg) —1 S)] SM n °t an appropriate reference value, because the coupling to the top 
quark (needed for the top quark triangle diagram driving the production through gluon-gluon fusion in the SM) is 
suppressed for extra scalars. As an example, let us consider the type t model of Eq. ([37]) . We have X_ = 0 and the 


( X + )tt coupling is suppressed by tp 1 for large tp. From Eqs. (26) we see that the first piece in the A^ H coupling is 
proportional to mtc a -p- The current LHC bounds on decays into vector bosons constrain s 2 a _p to lie within 20% of 
the SM value 1. Assuming the central value remains unity, Run 2 at the LHC should improve this bound to within 
5%. As a result, the gluon-gluon fusion production of H would be down with respect to the SM by roughly 0.05. 
Under these circumstances, the natural reference value for the gluon-gluon fusion production value would be 5% of 
the SM value. This can be accounted for by taking as the natural value r = 0.05. 


For a given value of r (say, 10%), Eq. ([62]) gives a limit on 


| a ma x |2 _|_ |&!” ax |-\ Qf coursej these values will be much 
larger than those obtained from Eq. (60). We would now like these values of |a^ ax | and |&^ ax | to survive the A ms, 
bound through the self-cancelation mechanism. We define the proportion of fine tuning required by 


Pft = 


bs 


-1K-/K, 


B s 


\ a bs \ 2 + l&, 


bs I 


7] Am 


b s 


ts / J., 

A a m Bs 


ms V 
fsj 


71 R (hs + Isb) 

E 2 ra [ pp(gg) -t S )] 


(63) 


SM 


where we have used Eqs. (§, (§, ©, and ( |62| ). For example, pft = 0.01 implies a fine tuning within 1%. A larger 
value of pft implies a lower fine tuning. Of course, the smaller the value chosen for r, the smaller the fine tuning 
(larger pft)- This is clearly seen in Fig. [4j showing pft as a function of ms for r = 1, r = 0.1, and r = 0.05. Consider, 
for example, a situation similar to that in the type t model, where the natural gluon-gluon production corresponds 
to r = 0.05, and take ms = 800GeV. We learn from Fig. [4] that for bs production to be equal to the gluon-gluon 
production requires a fine tuning in the self-cancelation within 1%. But this means that a reasonable self-cancelation 
within 10%, would imply a bs production to contribute around 10% of the total production. 

Notice that, in a generic model, the couplings of the spin 0 particle to bs , bd, and ds are independent. Thus, in a 
model independent search, each neutral meson system should be considered independently. However, in many models 
the couplings in each sector could have specific relations, as illustrated above in the BGL system. For definiteness, let 
us assume that the scalar S couples like one of the scalars in the type t model. According to Table [TTJ the couplings 
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Vs = 14 TeV 



FIG. 4: pft as a function of ms for r = 1 (solid line), r = 0.1 (dashed), and r = 0.05 (dotted). 


go like (X±)-.., which scale like the square of the entries of the matrix in Eq. (41). Namely, for the type t model 5 : 

\ /ji 1 r 


K 

m 2 s V 2 V 2 d *, 


B d 

m 2 b V 2 V t 2 d*, 


B s 

m 2 b V 2 V 2 *. 

(64) 


The key point is that both the scalar (a) and pseudoscalar ( b ) coefficients scale in the same way. From Eqs. 
we get 


Arns, 

A m Bd 


/I, m 'B B 

Ps K P+Hs K b\ 

PB d m B d 

a h K "‘ + 'fX( 


(65) 


Since K b /K 0 is almost the same in the Bd and B s systems and since, for each system, a and b scale in the same way, 
we expect that a 10% self-cancelation in the B s system will also imply a 10% self-cancelation in the Bd system. For 
the kaon system, besides the SM-like CKM rescaling, there is an m 2 s /m1 suppression. Indeed, taking 


C-sd 


m 2 s A 10 
A^ CBs 


( 66 ) 


for c = a,b, we find, 


A m K = ( — A 3 

,m b 


14 , K . 


K 


h ls K b] 


+ b i K ?° 


ffC m K 

Pb, m B s 


A nip 


(67) 


This means that a 10% self-cancelation in B s implies very small contributions to the mixing in the kaon system. 

Should it turn out that the self-cancelation mechanism is particularly effective (such that the FCNSI induced scalar 
production is a relevant percentage of the with gluon-gluon production), then one must look to other FCNSI effects 
for further constraints. As mentioned, decays of the type P° £ + £~ are only relevant if the scalar has a relevant 
coupling to t + 1~. In theories where this coupling is free, it could even vanish and no addition constraint arises. But in 
some theories this coupling is fixed by other quantities (such as masses and, possibly, elements of the leptonic mixing 
matrix), and they must be taken into account. For completeness, the appendix includes the relevant formulae. 

Finally, if the FCNSI couplings are large, they could induce S —» qiqj decays. Using the expressions in the appendix 
and the values of a max and & max in Table IV for ms = 125GeV, we find that the effects for S' —> bs + sb are less 


5 Notice that, in the type t mode, the CKM coefficients in Eqs. i |G4[ ) coincide with the coefficients that arise from the respective SM box 
diagrams. This means that the CKM multiplicative factor in the new contributions scale from system to system in accordance with 
experiment. 
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than a few percent of the total width, and completely irrelevant for the other FCNSI decays. As seen in Fig. [2] for a 
sufficiently large ms the maximum scalar partial width of S —> bs + sb could exceed S —> bb. But for masses above 
~ 135GeV, the decays into vector bosons take over. As a result, the S —> bs + sb decay will typically be a minute 
fraction of the total width. 


IV. CONCLUSIONS 

In general, a theory with more than one neutral scalar will induce FCNSI. Because these lead to mixing in the 
neutral meson systems, it is customary to eliminate or suppress such contributions. Strategies in the literature range 
from discrete symmetries - for example, a Zi symmetry in 2HDMs [37’, .38] -, to relation with CKM matrix elements 
[25] , and/or large scalar masses. In this article, we highlight a further possibility, which could occur for a spin 0 state 
with both scalar and pseudoscalar couplings. In that case, one could have a self-cancelation between both couplings. 

We explain in detail how this mechanism could occur even in models where the Higgs potential and the vacuum 
preserve CP. We illustrated this effect by showing the explicit couplings of the BGL models [23]. This case shows 
that a self-cancelation is quite natural, at least to some degree, depending on the exact values of the hadronic matrix 
elements. Our evaluations of this effect in the B systems are more promising using lattice estimates than using 
estimates with the vacuum insertion approximation. 

We also investigated the possibility that such FCNSI could induce a new mechanism of scalar production at LHC. 
In a theory with multi Higgs we expect that FCNSI can be ignored as a production mechanism, except in two cases: 
for some scalar where self-cancelation is active to some accuracy; or in limiting cases where the contribution from 
different scalars cancel each other, as illustrated in at the end of Section |HC| These effects should not be ruled out a 
priori. 
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Appendix A: Further FCNSI processes 

We are also interested in the decays of P° into two a pair of charged leptons P° —> £ + £~. We write the couplings 
of the scalar into the lepton pair as 


- C Y = S {£ (at + ibe 75 ) (-} -\ - 


(Al) 


where at and bt are real. We find 6 


T[P° -a i + r\ 


1 rn 5 P fl 

87 t (m q . + rn qi ) 2 nig 


\bji \ 2 Pe 


[l a t | 2 Pi + I bt\ 2 ] , 


(A2) 


where 


fo = 



(A3) 


and the form factor fp arises from 


(o\qj'y5Qz\P°) = i 


m. 


m q . + m qi 


fp, 


(0\q jqi \P°) = 0. 


(A4) 


Notice that a-, : does not appear in Eq. (|A2|) becaus e th e pseudoscalar meson P cannot couple to the scalar component 


of S, as seen in the second matrix element of Eq. (A4). 


6 The P° -4 P t: would seem to involve the different factor b. /;/ \. but, given Eq. 1M = M- 
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Finally, the Lagrangian in Eq. ([3]) induces h —» q^j decays. We find 


F [S -)• qiqj] 


N c 

8 n nig 


V a + a - [\ a ji\ 2a + + \ b Ji\ 2o! -] , 


(AS) 


where 


a± = m 2 s — (m 9i ± ) 


x 2 


(A6) 




